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Abstract 

5-H ■ . . 

We investigate the quantum Hall effect in graphene. We argue that in graphene in presence of an 
external magnetic field there is dynamical generation of mass by a rearrangement of the Dirac sea. 
We show that the mechanism breaks the lattice valley degeneracy only for the n = Landau levels 
and leads to the new observed = ±1 quantum Hall plateaus. We suggest that our result can be 
tested by means of numerical simulations of planar Quantum Electro Dynamics with dynamical 
fermions in an external magnetic fields on the lattice. 
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I. INTRODUCTION 



Graphene l|, |2|] is a flat monolayer of carbon atoms tightly packed in a two dimensional 
honeycomb lattice consisting of two interpenetrating triangular sublattices. The structure 
of graphene has attracted considerable attention from the fact that the low electronic exci- 
tations are described by the Dirac equation for two dimensional massless fermions {s, 4| (for 
recent reviews see Refs. 6|-[9| ). To describe the low-energy excitations one may consider the 
excitations at the Fermi level in the vicinity of two Dirac cones centered at two inequivalent 
corners of the Brillouin zone (called valley) K± = ± where Qq = 1.4210"^ cm IS 

the carbon-carbon distance. As a result, one obtains an effective Hamiltonian made of two 
copies of Pauli spinors which satisfy the massless two-dimensional Dirac equation with 
the speed of ligth replaced by the Fermi velocity vp — 1.0 10^ cm/s. In graphene, time 
reversal symmetry and parity, i.e., spatial inversion symmetry with respect to the center of 
a hexagon, transform electronic states between valleys. Thus, the energy spectrum is de- 
generate between the different valleys since the time reversal operation connects electronic 
states at to those at K_. If there is dynamical generation of gap Aq, then the energy 
levels of two valleys are related as ex+(Ao) = £:k_(— Aq) due to the inversion symmetry. 
In graphene immersed in transverse magnetic fleld, the relativistic massless dispersion of the 
electronic wave functions results in non-equidistant Landau levels (cgs units): 



2\n\H-^eH, n = 0, ±1, ±2 ... 



(I.l) 



where eH > 0, e being the elementary charge. The presence of anomalous Landau levels at 



zero energy, £0 = 0, leads to half- integer quantum Hall effect 
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where Rxy is the Hall resistance, Qs = 4: accounting for spin and valley degenerancy, and 
± stands for electrons and holes, respectively. This quantization condition corresponds to 
quantized fllling factor: 

1, 



iy = ±g,{N + 



±2, ±6, ±10 .. 



(1.3) 



Recently, experimental investigations of quantum Hall effect in graphene in very strong 
magnetic fleld H > 20T {IT = 10^ gauss) have revealed new quantum Hall states cor- 



responding to fllling factor u = , ±1 , ±4 



Id, 



111 ]. It is wiedely believed that the new 



plateaus at z/ = , ±4 can be explained by Zeeman spin splitting. On the other hand the 
z/ = ±1 plateaus are associated with the spontaneuos breaking of the valley symmetry in 
the n = Landau levels. Indeed, these states are naturally explained if there is dynamical 
generation of a gap Aq . Moreover, it turns out that Aq{H) ~ y/H 

This peculiar dependence on the magnetic field suggests that spontaneous generation of the 
gap is driven by the electron-electron Coulomb interaction jol [l2ll5|. In this picture, the 
dynamical gap is expected to be: 



Ao{H) ^ ^ ./^ ~ 163 ^KkB y/mT) (1.4) 



assuming the dielectric constant £ ~ 4 ll|]. The aim of this paper is to discuss a different 
mechanism which is able to spontaneuosly break the valley symmetry in presence of an 
external magnetic field. Indeed, we show that by a rearrangement of the Dirac sea the 
vacuum energy density is lowered by inducing a dynamical gap Aq ~ \/H. 
The plan of the paper is as follows. In Sect. Ullwe present the dynamical generation of a gap 
by a rearrangement of the Dirac sea in presence of an external magnetic field. In Sect. IIIII 
we discuss the fermion condensate in the one loop approximation. Finally, our conclusions 
are relegated in Sect. HVl 



II. DYNAMICAL GENERATION OF THE MASS GAP 

In this Section we discuss the mechanism which is able to spontaneously break the valley 
symmetry in presence of an external magnetic field by a rearrangement of the Dirac sea. In 
particular, we show that the vacuum energy density is lowered by inducing a dynamical gap 
Ao ~ VH. 

We said that the low-energy excitations in graphene are described by two Pauli spinors 
\I'-i- which satisfy the massless two-dimensional Dirac equation. When graphene is immersed 
in a transverse magnetic field, the relativistic massless dispersion of the electronic wave 
functions results in non-equidistant Landau levels Eq. (11.11) for both fields \E'-|- (see Fig. [H 
left). If there is dynamical generation of gap Aq, since the energy levels of two valleys 
are related as ex+iAo) = (— Aq) due to the inversion symmetry, the 'if± are no longer 
degenerate. Indeed, it is easy to find the following spectrum (see Fig. [H right): 

£+,o = Ao , e+,n = sign{n)^2\n\n^eH + A^ , n = ±l,±2..., (ILl) 
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FIG. 1: Schematic spectrum of Landau levels of graphene in applied magnetic field (left). Landau 
levels with dynamical generation of a gap Aq (right). The Fermi level is at ei? = 0. 



£_,o = -Ao = sign(n)y2|n|^^e// + Ag, n = ±1 , ±2 ... , (112) 

for the and fields respectively. 

Let us consider the vacuum energy density £ obtained by filling the Dirac sea. For massless 
fields, according to Eq. fll.ip we obtain: 



n=l 



where the factor in front of the sum takes care of the degeneracy of the Landau levels. In 
presence of a gap Aq, using Eqs. ( III. II) and ( III. 21) . we find: 



2nHc 
So that we get: 

AS = S{Ao) - S{0) 



n=l 



27r/ic 



-Ao -2 5^A/2n^%i^ + Ag +2 ^W2n^%/f ) . (n.5) 

n=l n=l 



-A + g(A) 




L A 



FIG. 2: Plot of the function —A + g{X) versus A. 



The sums in E(x_( III.5]) are divergent. To evaluate A£ we adopt the Schwinger's proper time 



regularization 



16|. To this end, suffice it to use the integral representation : 

ds d „„ 





a 



7CS ds 



A straightforward calculation gives: 



eH 



n^eH [ -A + giX) 



A 



Ar 



h^eH 



where: 



dx d 



TTX dx 



-2x 



1 - e 



-2x 



(11.6) 



(11.7) 



(11.8) 



In Fig. [2] we display the function —A + g{X)- It is evident that there is a negative 
minimum for A = A ~ 0.45. Thus, we see that it is energetically favourable to induce a 
dynamical gap: 

Ao =X \ fv^eH ~ 134 °K x ke ^H{T), 



(11.9) 



which is comparable to the estimate of the gap driven by the electron-electron Coulomb 
interaction, Eq. (II. 4p . 



III. CHIRAL CONDENSATE 



Recently, it has been obseved that the low-energy dynamics of graphene can be accounted 
for by Nf = 2 massless Dirac fields • These are descibed nonperturbatively by the lattice 
Euclidean action using staggered fermion fields [isl. As a consequence, we see that our 
mechanism for the dynamical generation of gap by a rearrangement of the Dirac sea can 
be tested by means of numerical simulations of planar Quantum Electro Dynamics with 
dynamical fermions in an external magnetic field on the lattice. In particular, the dynamical 
generation of the gap can be easily tested on the lattice by measuring the fermion condensate 
in the chiral limit. To see this, following Ref. [loj], we discuss relativistic fermions coupled 
with electromagnetic field in (2+1) dimensions. For reader convenience we do not follow the 
high-energy natural units where h = c = 1, but we employ cgs units. 

A spinorial representation in three dimensions is provided by two-component Dirac spinors. 
The fundamental representation of the Clifford algebra is given by 2 x 2 matrices which can 
be constructed from the Pauli matrices: 

7° = 0-3, 7^ = iai, 7^ = ia2 . (HI.l) 

Note that in this representation the mass term mc^^'^ is odd under both parity and time- 



reversal transformations (for a review, see Ref. 20|). We are interested in the chiral con- 
densate which is defined as: 



(^^) = 1- j (fx (0|^(a;)^(x)|0) , (III.2) 

where V is the spatial volume and |0) is the ground state (vacuum). In presence of an 
external constant magnetic field, it is easy to evaluate the chiral condensate in the one-loop 
approximation [l^: 



oo 



= _ mc^ > , III.3 

47r 27r ^ ^2nficeH + m?c^ 

It is worthwhile to stress that the first term in Eq. (1111.30 . due to the n = Landau levels, 
is independent on the mass and it is related to the parity anomaly. In fact, we see that the 
fermion condensate of the zero-mass theory in the presence of an external constant magnetic 
field corresponds to half-filled zero modes, which in turns leads to the observed half-integer 
quantum Hall effect in graphene. 



In three dimensions we have a reahzation of the Dirac algebra which is different from the 
2x2 representation Eq. fllll.ll) . In this representation the Dirac fermions are four-component 
spinor and the three 4x4 7- matrices can be taken to be j2l|: 

" ] y=f'"' ° ] .^-^i'"' " ] . (in.4) 

The representation given by Eq. (1III.4P corresponds to a reducible representation of the 
Dirac algebra. In this representation the fermionic mass term m^'$ is parity conserving. 
Indeed, comparing Eq. (IIII.ip with Eq. (IIII.4P we see that a four-component fermion with 
mass m corresponds to two two-component fermions with mass m and — m respectively. So 

that, if we write \1/ = I ' I , we readily get: 



(^^) = (^+$+ - ^_^_) . (III.5) 



Since under parity transformation ^+ — t- and — t- \1'+, we see that the chiral conden- 
sate Eq. (jIII.51) is parity invariant. In presence of an external constant magnetic field, we 
may easily evaluate the chiral condensate in the one-loop approximation by using Eqs. fllll.SP 
and f lTTLSll : 



(^^) = -2|m|c^^y , ^ . (III.6) 

27r ^ ^2nficeH + m?c^ 



Using the following identity: 

1 r ds 



we recast Eq. (1III.6P into: 



e""' , (III.7) 



-ax 



,— , , r, HceH 1 dx „,2^ e 

(\1/\1>) = -2 mc , / ^= e 

27r ^2nceH Jo 



a = ' . (III.8) 

yj2hceH 



Now, we may regularize the integral by analytic continuation. To do this, we note that |22| : 







v-l -^ix 1 

^-13^ = ^r(.)c(.,^), (III.9) 



where T{x) is the Euler gamma function and ({x,y) is the generalized Riemann Zeta func- 
tion. Combining Eqs. f llll.SP and fllll.Qp we finally obtain: 



HceH 1 r(i) ^ 1 



,2 



( W) ^-2Hc^——^a.,,,a^U^^-^. („I,0) 
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Our results can be applied to graphene if the Coulomb interaction can be neglected. In that 
case, observing that Aq = |m|c^ we readly obtain from Eq. flIII.10|) : 

/- ^ . ficeH 1 r(i) ^1 . 



a = . (III.ll) 

2H^eH 

For small gap, we may expand to the first order in Aq. Using r(i) = ^/^T and ({x, 1) = C{^)y 
we get: 

fire H 1 1 
m = -2Ao^^i=4). (III.12) 

This last equation can be used to extract the gap Aq from non perturbative measurements 



of the fermion condensate on the lattice. Indeed, preliminary Monte Carlo simulations [23 1 
of planar Quantum Electro Dynamics with two degenerate staggered fermions on the lattice 
suggest that external magnetic fields give rise to the spontaneous breaking of the chiral 
symmetry in the continuum limit. 



IV. CONCLUSIONS 

We have discussed the spontaneous generation of the gap in graphene induced by a 
rearrangement of the Dirac sea. We find that the gap Ao(H) ~ \/H, and it is comparable to 
the gap driven by the electron-electron Coulomb interaction. We argued that our dynamical 
generation of the gap can be easily tested on the lattice by measuring the fermion condensate 
in the chiral limit. 



[1] K. S. Novoselov et al, Science 306 (2004) 666. 

[2] K. S. Novoselov et al, Proc. Natl. Acad. Sci. U.S.A. 102 (2005) 10451. 
[3] K. S. Novoselov, A. K. Geim, et al, Nature 438 (2005) 197. 
[4] Y. Zhang et al, Nature 438 (2005) 201. 

[5] K. S. Novoselov, E. McCann, et al, Nature Physics 2 (2006) 177. 
[6] A. K. Geim and K. S. Novoselov, Nature Materials 6 (2007) 183. 
[7] M. I. Katsnelson, Materials Today 10 (2007) 20. 

8 



[8] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. 
Phys. 81 (2009) 109. 

[9] V. N. Kotov, B. Uchoa, V. M. Pereira, A. H. Castro Neto, and F. Guinea. Electron- Electron 
Interactions in Graphene: Current Status and Perspectives, (2010) I ar Xiv :1012.3484l 
[10] Y. Zhang, et al, Phys. Rev. Lett. 96 (2006) 136806. 

[11] Z. Jiang, Y. Zhang, H. L. Stormer, and P. Kim , Phys. Rev. Lett. 99 (2007) 106802. 
[12] D. V. Khveshchenko,Phys. Rev. Lett. 87 (2001) 206401. 

[13] V. P. Gusynin, V. A. Miransky, S. G. Sharapov, and 1. A. Shovkovy, Phys. Rev. B 74 (2006) 
195429. 

[14] L F. Herbut, Phys. Rev. B 75 (2007) 165411. 

[15] V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, Int. J. Mod. Phys. B 21 (2007) 4611. 
[16] P. Cea, Phys. Rev. D 32 (1985) 2785; Phys. Rev D 34 (1986) 3229. 

[17] J. E. Drut and T. A. Lahde, Phys. Rev. Lett. 102 (2009) 026802; Phys. Rev. B 79 (2009) 
165425. 

[18] S.J. Hands, J.B. Kogut, C.G. Strouthos, Nucl. Phys. B 645 (2002) 321. 
[19] P. Cea and L. Tedesco, J. Phys, G 26 (2000) 411. 

[20] R. Jackiw, Massive gauge theories in three- dimensions (= at high temperature) presented at 

Artie School of Phys., Akaslompolo, Finland, 1982. 
[21] T. Appelquist, M. Bowick, D. Karabah, and L.C.R. Wijewardhana, Phys. Rev. D 33 (1986) 

3704. 

[22] LS. Gradshteyn and LM. Ryzshik, Table of Integrals, Series and Products, 1980 (New York, 
Academic) . 

[23] P. Cea, L. Cosmai, P. Giudice, and A. Papa, PoS LATTICE2011 307 (2011). 



9 



